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Novel imaging spin-filter techniques, which are based on low energy electron diffraction, are cur¬ 
rently of high scientific interest. To improve the spin-detection efficiency a variety of new materials 
have been introduced in recent years. A new class of promising spin-filter materials are represented 
by multiferroic systems, as both magnetic and electric ordering exist in these materials. We have 
investigated Fe/BaTiO3(001), which defines a prominent candidate due to its moderate spontaneous 
polarization, for spin filter applications calculating diffraction patterns for spin polarized electrons 
incident on the Fe surface. Motivated by the fact that spin polarized low energy electron diffraction 
is a powerful method for the determination of the properties of surfaces we investigated the influ¬ 
ence of switching the BaTiOs polarization on the exchange and spin orbit scattering as well as on 
reflectivity and figure of merit. This system obviously offers the possibility to realize a multiferroic 
spin filter and manipulating the spin-orbit and exchange scattering by an external electric field. The 
calculations have been done for a large range of kinetic energies and polar angles of the diffracted 
electrons considering different numbers of Fe monolayers. 


I. INTRODUCTION 

Materials combining different ferroic properties are 
called multiferroics [l|. Multiferroic systems, especially 
multiferroic heterostructures, are promising in techni¬ 
cal applications giving for example the opportunity to 
change the magnetization of a ferromagnetic material al¬ 
tering the electric polarization of a ferroelectric material 
or vice versa. An essential requirement of the realiza¬ 
tion of such systems was the development of experimen¬ 
tal techniques over the last years for growing epitaxial 
magnetic and ferroelectric materials Hi- Additionally 
a well-defined interface has to be ensured between both 
ferroic phases to enable an effective coupling. Various 
experimental groups investigated prototypical devices us¬ 
ing external electric fields for a change of the ferroelectric 
polarization and investigated the influence on the mag¬ 
netic properties of the ferromagnetic phase Q. The ap¬ 
plication of such new systems are broad but especially as 
memory devices epitaxial thin films might be interesting 
(Ef. Because of the technical applicability the discovery 
for new multiferroic materials with strong coupling be¬ 
tween the ferrom^netic and the ferroelectric phases has 
been strengthen [El, 0- The understanding of the cou¬ 
pling mechanism between both ferroic phases is essential 
for the development of multiferroic heterostructures and 
their technical application. Therefore different ferroelec¬ 
tric materials have been studied extensively in previous 
theoretical investigations where the BaTiOs (BTO) takes 
a prominent position due to its moderate spontaneous 
polarization [8[. 

For the construction of a multi-channel vector spin po- 
larimeter at BESSY II new types of material classes are 
tested for their applicability as reflecting mirror for spin 
polarized electrons. Beside well tested materials likewise 
W, Ir, and oxygen passivated Ee other material classes 


should be investigated yielding interesting properties for 
spin filter purposes. Because of the high magnetic mo¬ 
ments for the Ee layers of Ee/BTO and the d^endency 
of their magnitude on the BTO polarization [y, this 
heterogeneous system is an interesting candidate for the 
application as spin polarizing mirror. In principle using 
a multiferroic heterostructure would give the possibility 
to alter the exchange and spin-orbit scattering switching 
the electrical polarization of the BTO. Thereby no addi¬ 
tional magnetic field is necessary for changing the surface 
magnetization of the ferromagnetic material. Especially 
working areas inducing an altering of the scattered spin 
direction by switching the BTO polarization are interest¬ 
ing for an application as spin polarizing mirror. 

In this work we studied the scattering of spin polarized 
electrons from the surface of the multiferroic heterostruc¬ 
ture Ee/BTO. We calculated diffraction patterns visual¬ 
izing the dependency of the scattered intensity on the 
polar angles and kinetic energies. Because the polariza¬ 
tion of the electrons were oriented perpendicular to the 
scatte ring pla ne, both exchange and spin-orbit scattering 
occur [nl, [ 13 . Therefore we investigated the so-called ex¬ 
change and spin-orbit asymmetries as well as the effective 
reflectivity and the figure of merit (EOM). Additionally 
we investigated the layer dependence of the exchange and 
spin-orbit scattering. We applied SPLEED (Spin Polar¬ 
ized Low Energy Electron Diffraction) calculations to 1, 
2 and 3 monolayer (ML) Ee on top of BTO because their 
applicability as multiferroic system has been mentioned 
and their electronic and magnetic structure was investi¬ 
gated in detail in previous works Bill. Eor all systems 
the crystal structure was taken from previous investiga¬ 
tions providing a relaxed interface and surface [9[. This 
accounts for the fact that SPLEED is very surface sensi¬ 
tive due to the low kinetic energies of the incident elec¬ 
trons. Eor the calculation of the electronic properties and 
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the electron diffraction we applied the fully relativistic 
KKR-method in the framework of spin-polarized density 
functional theory to account for effects based on exchange 
and spin-orbit interactions in one step Hi Hi- 

The paper is organized as follows: In Sec. [II] we de¬ 
scribe the application of the underlying SPLEED theory 
in the framework of the SPR-KKR method. In Sec. IIIII 
we discuss our computational results concerning the elec¬ 
tron diffraction. In Sec. IE] we summarize our results. 


II. THEORETICAL APPLICATION 

We briefly introduce the theoretical method imple¬ 
mented in the SPR-KKR program package for the calcu¬ 
lation of spin-polarized electron scattering from arbitrary 
surface systems. The SPLEED calculations are done us¬ 
ing the layered-KKR method 0- This method describes 
the scattering of spin polarized electrons from a stack of 
atomic layers representing a semi-infinite surface system. 
The two main steps are the treatment of multiple scat¬ 
tering within one specific atomic layer and the scattering 
between the atomic layers of the stack. The combination 
of both scattering mechanisms result in the calculation 
of the so-called bulk reflection matrix. With the deter¬ 
mination of the bulk reflection matrix the diffraction of 
spin-polarized electrons from a surface is defined 00. 
Considering exchange and spin-orbit scattering there ex¬ 
ist mainly two configurations according to the orientation 
of magnetization and polarization of the electron with 
respect to the scattering plane. Whereby the scattering 
plane is spanned by the wave vector of the incident elec¬ 
tron k and the scattered electron k' (see Eig. [T]). In our 
calculations both surface magnetization and polarization 
of the electron are oriented perpendicular to the scatter¬ 
ing plane. Therefore exchange and spin-orbit scattering 
contribute 0. 

The calculation of the scattered electron intensity was 
done for all combinations of magnetization and polar- 
ization, i.e. one ends up with four different intensities 
) [ill, uM • Here /i is the magnetization direction and 
a gives the polarization direction of the electron. Erom 
the reflected intensities the exchange asymmetry, spin- 
orbit asymmetry and EOM are calculated. The spin-orbit 
asymmetry can be determined via the relation [hI 

Aoc = 1{A+-A.), ( 1 ) 

which can be ascribed to the reflected intensities using 
the definitions 
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The exchange scattering for a specified magnetization 
direction is defined as A+ and A_. 

The most important quantity for the characterization 
of the diffraction of electrons from a surface is the EOM. 
The reflectivity as well as the asymmetry (exchange or 
spin-orbit) contribute to this quantity via the following 
equation 


FOM+(_) =/+(_). 4(_). (4) 

In Eq. (H]) the EOM is calculated separately for the 
two specific magnetization directions. This gives insight 
into changes of the scattering behaviour changing the 
magnetization direction at the surface. The spin orbit 
induced EOM is defined by 


■ Al„ (5) 

with leff the effective reflectivity 


Ieff=lt + I++ll+IZ. (6) 

Eor the calculation of SPLEED patterns it is nec¬ 
essary to determine the single site scattering matrix 
for the individual atomic types involved in the half¬ 
space of Ee/BTO. The self-consistent potentials neces¬ 
sary for calculating the single-site scattering matrices 
were taken from previous works @,0. With the sin¬ 
gle site-scattering matrix we were able to calculate the 
Kambe X-matrix and the bulk reflection matrix as de¬ 
scribed above. 

Considering the escape of the diffracted electrons into 
vacuum an essential quantity for SPLEED is the work 
function and the surface potential barrier. Eor the work 
function we applied 4.7 eV which is a reasonable value 
for Ee(OOl) [l3 whereas the surface potential barrier was 
simulated by the Rundgren-Malmstrom parametrization 
0. We calculated SPLEED patterns for a broad range 
of kinetic energies and polar angles according to the 
working areas as scattering mirror. The calculations were 
done for the specular reflected beam and for the polariza¬ 
tion of the BTO pointing in the direction of the surface 
normal (Pnp) and to the opposite direction (Pdown)- 

In Eig. [Tjthe side view and top view of the Ee/BTO 
heterostructure are shown. Here 3 ML Ee are placed on 
top of the substrate BTO. The unit cell of BTO corre¬ 
sponds to a tetragonal distorted structure (P4mm) with 
a lattice parameter of 3.943 A0. The structural details 
have been published in previous works i0- 
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FIG. 1: Half-space of 3 ML Fe/BTO used for the SPLEED 
calculations. Left: Side view. Right: Top view. The vari¬ 
ous atomic types have been indicated. The reflected electron 
intensity (7^(0, $,E)) is drawn. 0 is the polar angle of the 
incident electrons and $ the azimuthal angle. Additionally 
the parameters which have been varied are shown: The elec¬ 
tric polarization of BTO (Pup, ^down), the polarization of 
the electron (a), the in plane magnetization of the surface 
(/i). The plane of incidence is spanned by the wave vectors 
{k,k') and is aligned along the [010]-direction. The polariza¬ 
tion of the electron and the surface magnetization are altered 
along the [100]-direction. 


Remarkably the surface of the BTO has a (001) orien¬ 
tation and is terminated by O and Ti. The termination of 
the Fe/BTO system has been investigated in previous ex¬ 
perimental studies [ 2 ^ . The multiferroic heterostructure 
was simulated by setting up a half-space of four unit cells 
of BTO for the transition from the surface to the bulk. 
From the fifth unit cell on the bulk region was started. 
Due to the distorted unit cell of BTO (P4mmm) a re¬ 
manent electrical polarization occurs. This result from a 
shift of the Ti atom and the O atoms in opposite direction 
and generates the ferroelectric properties of BTO [ 3 . It 
has been shown that the first Fe layer (Feo) on top of 
BTO is positioned on top of the O atoms. The second Fe 
layer has two inequivalent Fe sites which are on top of Ti 
(Fexi) and Ba (FcBa)- The third Fe layer (Fej) is placed 
on top of the Fe atoms of the first Fe layer [9|. The 
surface as well as the interface of the multiferroic het¬ 
erostructure has been relaxed in previous works H u sing 
the Vienna Ab initio Simulation Package (VASP) 

The resulting structure parameters have been presented 
and discussed Using the relaxed crystal structure 
we applied the fully relativistic multiple scattering for¬ 
malism for the calculation of the SPLEED patterns as 
sketched above M- The self-consistent calculation of 
the system potentials has been done using a previously 
multicode approach [13, IH, . 


III. RESULTS AND DISCUSSION 

The results of the electronic structure calculations, es¬ 
pecially the spin magnetic moments have shown to be in 
good agreement with published data i,®. The mag¬ 
netic properties of the Fe surface are most important 
for the exchange scattering of spin-polarized electrons. 
Changing the polarization of the BTO affects the mag¬ 
netic moments of the Fe layers. Therefore the surface 
magnetic moments have shown to react on a competition 
of the reduced number of nearest neighbors at the sur¬ 
face and hybridization effects between the Fe 3d and O 2p 
states as well as the Fe 3d and Ti 3d states [H, [ol, [25| . The 
impact on the electronic structure has been investigated 
in detailed studies [9[ . It was shown that for 1 ML Fe on 
top of BTO a ferromagnetic ground state occurs with a 
large spin magnetic moment for Fe due to the reduced co¬ 
ordination number at the surface. For 2 ML Fe on BTO 
a ferrimagnetic ground state was predicted whereas the 
Fe atoms in the second Fe layer (Fexi, Feea) show an an¬ 
tiparallel alignment of their spin directions and different 
magnitudes. For 3 ML Fe a ferromagnetic ground state 
was predicted. The change of the BTO polarization af¬ 
fects the spin magnetic moments of Fe via a change in the 
hybridization of the Fe, Ti and O states mentioned above. 
A smaller (larger) hybridization of the 3d states of Fe - 
Ti and the 2p states of O results in a larger (smaller) spin 
magnetic moment of Fe 0 - Therefore the magnetoelec¬ 
tric coupling between the ferroelectric and ferromagnetic 
material is based on hybridization effects [9[ . The change 
of the magnetic properties has an influence on the opti¬ 
cal properties like the absorptive part of the conductivity 
tensor [ 1 ^ and are detectable using the scattering of elec¬ 
trons with low kinetic energy, as will be shown below. 

In Fig. [ 2 ] the reflectivity for 3 ML Fe/BTO for a [100] 
surface magnetization is shown. 



Ekin (eV) (eV) 


FIG. 2: Diffraction pattern of the reflectivity for 3 ML 
Fe/BTO for a [100]-orientation of the magnetization direction 
(color online). Left: The polarization of the BTO is directed 
along the surface normal (P^p)- Right: The polarization of 
the BTO is directed in opposite direction (Pdown)- 

The range of polar angles and kinetic energies have 
been chosen according to a possible application as scat¬ 
tering mirror [26|. The diffraction pattern can be di- 
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vided into two main parts with respect to the kinetic en¬ 
ergy. The main differences occur for kinetic energies < 8 
eV, i.e. below the emergence threshold. The emergence 
threshold marks the occurrence of an additional beam 
lowering the intensity of the specular diffracted one. The 
area of high reffectivity for kinetic energies at 4 eV covers 
nearly the total range of polar angles. This is important 
for later application as scattering mirror giving the pos¬ 
sibility to vary the position of the mirror with respect 
to the polar angle of the incident electrons. For kinetic 
energies lower the emergence threshold a large sensitivity 
changing the BTO polarization is visible. On the other 
hand for kinetic energies above the emergence threshold 
switching the BTO polarization has less inffuence on the 
reffectivity. This is due to the fact that electrons with 
lower kinetic energy react more sensitive on changes of 
the surface magnetization, i.e. onto variations of the ex¬ 
change and spin-orbit scattering. For kinetic energies 
above 10 eV an area of high reffectivity is visible for po¬ 
lar angles around 30°. Because of the high kinetic energy 
this area is much less affected by a change of the BTO 
polarization. 

In Fig. [3] the exchange asymmetries for 3 ML Fe/BTO 
for both BTO polarizations are shown with an [100] 
magnetization direction of the surface. For the investi¬ 
gated range of polar angles and kinetic energies two areas 
(marked by black rectangles) occur with a large variation 
of the exchange scattering. 

The first area is located at lower kinetic energies 
of the incident electrons (< 6 eV) the second area at 
higher kinetic energies (12-17 eV). The different colors 
in the diffraction patterns indicate a change of the align¬ 
ment of electron polarization and surface magnetization. 
Whereas the red color correspond to a parallel alignment 
of polarization and magnetization the blue color corre¬ 
spond to an antiparallel alignment. As shown in Fig. 
[3] changing the BTO polarization from Fup to Fdown 
changes the refiected spin polarization in the marked 
areas. Especially for the area at lower kinetic energy 
and higher polar angle, which provides high reffectivity, 
a change of the refiected spin direction due to a change 
of the BTO polarization is possible. 

In Fig. [3] the exchange scattering for both polariza¬ 
tion directions of the BTO for a magnetization in [100] 
direction is shown in addition. In comparison to a mag¬ 
netization direction pointing in the [100] direction a less 
pronounced change in the exchange scattering occurs es¬ 
pecially for higher kinetic energies. This results from 
spin-orbit scattering caused by the fact that for a scat¬ 
tering configuration with vanishing spin-orbit interaction 
A+ = —A- holds Ha, i.e. the exchange asymmetries for 
both magnetization directions are indistinguishable. 

The spin-orbit induced contribution to the specular 
scattering of the electrons is shown in Fig. IH Marked by 
black rectangles the spin-orbit asymmetry shows a higher 
change by switching the BTO polarization. These areas 
coincide with the rectangles for the exchange asymme¬ 
try for a [100] direction of the magnetization. The blue 
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FIG. 3: top panel: Diffraction pattern of the exchange 
asymmetry (A+) for 3 ML Fe/BTO. The magnetization of 
the surface is along the [100] direction (color online). Left: 
The polarization of the BTO is directed along the surface 
normal (Pnp). Right: The polarization of the BTO is directed 
in opposite direction (Pdown)- bottom panel: Diffraction 
pattern of the exchange asymmetry (A-) for 3 ML Fe/BTO. 
The magnetization of the surface is along the [TOO] direction 
(color online). Left: The polarization of the BTO is directed 
along the surface normal (Pnp)- Right: The polarization of 
the BTO is directed in opposite direction (Pdown)- 


rectangles correspond to an exchange asymmetry for the 
[100] direction. According to the specified areas of the 
exchange asymmetry the blue rectangles show a minor 
change of the spin-orbit asymmetry and therefore a mi¬ 
nor change in the exchange asymmetry (see Fig. [3j). This 
results in a more sensitive reaction of the exchange asym¬ 
metry on a change of the BTO polarization using a [100] 
orientation of the magnetization. This is important for 
the use as spin polarizing mirror indicating that an orien¬ 
tation of the surface magnetization would be preferable 
along the [100] direction. 

The lowering of the reffectivity affects the FOM via 
Eq. dH). The exchange asymmetry enters Eq. (j4]) to the 
power of two and therefore dominantly determines the 
EOM. This is the reason for the correspondence of areas 
of high exchange scattering and EOM. In Eig. [5]the EOM 
for both BTO polarization directions are shown. 

The EOM shows significant changes for kinetic ener¬ 
gies < 6 eV. Therefore the multiferroic heterostructure 
3 ML Ee/BTO would have most promising technical ap¬ 
plicability for VSPLEED experiments. With respect to 
the area of high EOM around 4 eV for an polarization of 
the BTO oriented along the surface normal the system 
Ee/BTO could be used as switch for the electron spin 
directions. Changing the BTO polarization the value of 
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FIG. 4: Diffraction pattern of the spin-orbit asymmetry 
(Asoc) for 3 ML Fe/BTO (color online). Left: The polariza¬ 
tion of the BTO is directed along the surface normal (Pnp)- 
Right: The polarization of the BTO is directed in opposite 
direction (Pdown)- 
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FIG. 5: Diffraction pattern of the FOM for a [100] orienta¬ 
tion of the magnetization direction for 3 ML Fe/BTO (color 
online). Left: The polarization of the BTO is directed along 
the surface normal (Pup)- Right: The polarization of the 
BTO is directed in opposite direction (Pdown)- 


the FOM drops drastically, i.e. no electrons are reflected. 

In comparison to investigations done for an oxygen pas¬ 
sivated Fe surface HE 113 the reflectivity is in the same 
order whereas the FOM is an order of magnitude lower 
for the multiferroic heterostructure. Therefore the spin 
resolving power is less for the Fe/BTO system lowering 
its applicability as spin-polarizing mirror. Because it is 
know that Fe(001)-p(lxl)-0 has a longer lifetime in vac¬ 
uum and additionally a higher FOM in comparison 
to a Fe(OOl) surface it would be interesting to investigate 
the effect of an O overlayer on top of Fe/BTO which will 
be investigated in a later work. 


A. Layer dependence of SPLEED 


of the BTO occur for 3 ML Fe/BTO. Therefore the re¬ 
sults for 1 ML and 2 ML will not be shown in detail 
here. Nevertheless, the scattering patterns for 1 ML 
and 2 ML reveal differences between both BTO polar¬ 
izations corresponding to 3 ML Fe/BTO. The differences 
in the exchange scattering between 1 ML, 2 ML and 3 
ML Fe/BTO can be related to the changes of the in plane 
spin magnetic moments changing the polarization of the 
BTO. In Table [I] the changes of the in plane spin mag¬ 
netic moments for a change of the BTO polarization are 
listed for the several Fe layers. The largest change of the 
spin moments occur for the topmost Fe layer of 3 ML 
Fe/BTO. Therefore an electron impinging on the sur¬ 
face interacting mainly with the topmost Fe layer would 
adept the largest change of the exchange scattering for 3 
ML Fe/BTO. For 1 ML and 2 ML Fe/BTO the changes 
are smaller resulting in a less pronounced change of the 
exchange scattering. 



1 ML 2 ML 3 ML 

Eei 

-0.18 

EeTi+Eesa 

-0.08 0.04 

Feo 

-0.02 0.0 0.2 


TABLE I: Layer resolved change of the in plane spin magnetic 
moments (Arris = rns(Pup) - rns(Pdown)) by changing the 
polarization of BTO. Units are given in [fis]- 


IV. SUMMARY 

Using an ab initio formalism we calculated SPLEED 
patterns for 1, 2 and 3 ML Fe/BTO. We investigated 
the change of the reflectivity, the FOM and both ex¬ 
change and spin-orbit asymmetry changing the polariza¬ 
tion of BTO. We showed that the largest differences of the 
diffraction patterns by changing the electric polarization 
of BTO results for 3 ML Fe/BTO. We located two areas 
in the diffraction patterns for which a change of the BTO 
polarization induces a significant change of the exchange 
asymmetry. Therefore a combination of a ferroelectric 
and a ferromagnetic material would enable new applica¬ 
tions for spin filter purposes, according to the working 
areas especially for VSPLEED experiments. For higher 
kinetic energies the Fe/BTO system has a smaller FOM 
compared to the intensively studied Fe(001)-p(lxl)-0 
surface. Therefore investigating the impact of an ad¬ 
ditional O layer on top of Fe will be done in the near 
future. 


In addition we investigated the change of the electron 
scattering for thinner Fe layers (1 ML, 2 ML) on top of 
BTO. The structural properties correspond to systems 
studied in previous works [9|. The results of our calcu¬ 
lations have shown that the most pronounced changes 
of the exchange asymmetry by changing the polarization 
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